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Abstract 

We investigate the finite size corrections to the equilibrium magnetization of an 
Ising model on a random graph with N nodes and N"' edges, with 1 < 7 < 2. 
By conveniently rescaling the coupling constant, the free energy is made extensive. 
As expected, the system displays a phase transition of the mean-field type for all 
the considered values of 7 at the transition temperature of the fully connected 
Curie- Weiss model. Finite size corrections are investigated for different values of the 
parameter 7, using two different approaches: a replica-based finite expansion, and 
a cavity method. Numerical simulations are compared with theoretical predictions. 
The cavity based analysis is shown to agree better with numerics. 
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1 Introduction 



Complex networks of interacting elements are ubiquitous in nature and display 
a rich and fascinating phenomenology which remains to be fully elucidated 
pH2] . Coherent macroscopic patterns can emerge from a limited set of rules 
which govern the microscopic evolution of the elementary constituents. Besides 
the specificity of the assigned interactions {e.g. long versus short range), global 
coherence reflects also the peculiarity of the underlying network of connections. 
An emblematic, though extreme, example, is the fully coupled network, where 
each one of the N nodes experiences a direct link with all the others. In the 
opposite limit, the number of links per node is independent of the size of the 
network. Studying the modifications induced on the statistical properties of 
the system by such dilution has important implications for a broad class of 
applications. 

Most studies have so far focused on the finite connectivity scaling, that is 
graphs where the number of nodes and the number of links Nl go to infinity, 
keeping the ratio N/Ni fixed. Models defined on such graphs have been exten- 
sively studied in the physics literature, especially in the context of spin glasses 
and combinatorial optimization [3|1] . Several recent mathematical papers also 
analyze rigorously ferromagnetic Ising models on such graphs [5|6|7ll8]|9] . 

Besides the fully coupled case, comparatively very few papers have been de- 
voted to the infinite connectivity scaling, where A^, A^^, and Nl/N all go to 
infinity. Bovier and Gayrard [10] showed rigorously that, under an appropri- 
ate rescaling, the Ising model defined on such graphs is completely equivalent 
in the infinite A^ limit to its fully connected [Nl = N{N — l)/2) counter- 
part. On the basis of this finding, it can be safely conjectured that analogous 
conclusions should apply to other spin models when inspected on a similar 
geometry. 

The rigorous result by Bovier and Gayrard leaves however open the question of 
the speed of convergence towards the fully connected limit. Finite size effects 
are indeed crucial when one aims at understanding the behavior of a finite 
graph. Given a finite random graph with relatively high connectivity, two 
different approaches may be used to address this issue: on one hand, one might 
consider expanding around the fully connected solution of ^0], to compute the 
leading order finite A^ corrections; on the other hand, the graph under scrutiny 
can be seen as a realization of an Erdos-Renyi random graph [TT], with high. 
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but finite, mean connectivity. This observation enables one to employ the 
powerful techniques developed for this latter scaling. 

In this paper, we will perform both computations. The first by resorting to 
the celebrated replica trick [12]; the second, by means of the cavity method 
|13] . Theoretical estimates will be then compared with numerical simulations, 
performed for different sizes and connectivities, in order to test their accuracy. 

The paper is organized as follows. In Section [21 we introduce the model and 
present the infinite solution. Numerical simulations are presented in Sec- 
tion [3l where the role of finite size corrections is quantified, with reference to 
a selected, macroscopic observable (magnetization). Section [His dedicated to 
developing the finite size replica calculations, whose predictions are compared 
with simulations. The subsequent Section [5] presents the cavity based analysis. 
Finally, in Section [H] we sum up and conclude. Two Appendices, devoted to 
the annealed model and to the Curie- Weiss model, are added in order to give 
the technical details on the textbook developments of Section [2l 



2 The model 



We here consider an Ising model defined on a uniform random network topol- 
ogy. The network is made of sites, each tagged by a discrete counter i, which 
ranges from 1 to N. On each site sits an Ising spin variable S'j = ±1. Two ran- 
domly selected nodes, say i and j, are connected through a coupling constant 
Jij. The number of links A'^^ is bounded from above by = A^(A^ — l)/2, 
which corresponds to the fully connected case, since we are avoiding double 
edging of two sites and self wiring. We introduce the dilution parameter 7 by 
scaling the number of links as = 

= iVT/7! (1 + 0(1/A^)), the normal- 
ization factor 7! is introduced so that the fully connected topology is exactly 
reproduced when 7 — > 2. We restrict our analysis to the interval 1 < 7 < 2, 
which corresponds to network topologies that range from a number of links 
growing linearly with the size of the system (7 = 1) up to the fully connected 
case (7 = 2). 

The Hamiltonian of the Ising model on such a diluted network is 




With this scaling of the coupling constant the energy is extensive. In the 
simplest formulation, Jij is set to an identical reference value, say J > 0, if 
the nodes i and j are connected, zero otherwise. In the following we shall 
study the behaviour of the system as a function of the dilution rate 7. As 
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proved in [TO], for 7 strictly larger than its lower bound 1, a second order 
phase transition of the Curie- Weiss type always occurs. Our main goal is to 
characterize the behaviour of the system for a finite size A^. In the following 
we start by analyzing the large N limit. The partition function Z reads 



Z = ^ exp JijSiSj) = E n exp (-^■J^,3S^SJ] 

{5J ^^^Li^j {S,}i^j ^^^L 



(2) 



where /? stands for l/ksT, where ks is the Boltzmann constant and T the 
system's temperature. The outer sum in Eq. ([2]) extends over all possible spin 
configurations. The coupling factor Jij is related to the linking probability via 
the following condition 



J with probability p 
Jij = \ (3) 
with probability 1 — p , 



where p is 



p = N,/N = ^N^'-^ (1 + (1) ) . (4) 



The 0{1/N) term in the probability represents subdominant finite size effects, 
and plays no role in the following. The probability distribution P{Jij) reads 

P(J,^)=p6iJ,,-J) + il-p)6iJ.,) . (5) 

We end this section by showing that for all 1 < 7 < 2, this system is exactly 
equivalent, in the N —>■ 00 limit, to the fully coupled Curie- Weiss model. This 
is a non rigorous rephrasing of the main result in Ref. [10], which will set the 
stage for the finite- studies of the following Sections. 

We want to compute (InZ)j, where denotes the average over disorder. 
This is achieved via the celebrated replica trick, which is based on the identity 

(lnZ), = lim J^"^"~^ , (6) 

where n is a assumed to be a real number. The central idea consists in carrying 
out the computation for all integers n, extending the results for all n, and 
performing in the end the limit for n — > 0. 

In our setting the replicated partition function reads 
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Z"=[^exp(-/?i7) 

{s-} „ 



where the index a runs over the n rephcas. 



Averaging over the disorder returns 



{z-)j=j:p{j.j) e exp I tf^^E-m 



{Sf} 



2'' m-^ f-"''" 



{Sf} Jij 



where Tij = J^a^iSj. RecaUing Eq. ([5]), one straightforwardly obtains 



in J = E n [i - --^ + 4^e^^i^^'^ 



7! iV2-7 -y! iV2-7 



To proceed further, we now expand the exponential function to its 
order approximation, which immediately yields 
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and, expanding the logarithm 



)i = E (^E + - 2m^^^ 



Keeping only the leading order in N, we have 




where Zmf is the partition function of the Curie- Weiss model. Thus, we re- 
cover the fact that at leading order in A^, the dilute model is equivalent to the 
fully coupled one for all 1 < 7 < 2. The calculation of Zmf is standard and 
recalled in Appendix B. We summarize here the results for convenience. The 
magnetization of the system m = limTv^oo J2iSi/N is obtained by solving the 
implicit equation 

m = tanh(2/3m) . (16) 

The critical inverse temperature, jSc = 1/2, separates the non magnetized 
phase from the magnetized one. The probability distribution function of the 
magnetization can be calculated from the free energy F{m) as 

^^"^^ = V(W\ ^ (-/3F(m)Ar) , (17) 

^yPl {C\m{C)=m} 

where C represents the subsets of spins configurations that have magnetization 
equal to m. For small m, one can expand F{m) in powers of m and obtain 

P{m) oc exp (^—€2^12^ — C4m^^ . (18) 

At the critical point (3c = V^^ one gets C2 = and C4 = 1/12. We stress once 
again that all these results do not depend on 7 in the N —>■ 00 limit. 

In the following Section, we will discuss the numerical implementation, test 
the above infinite A^ theory and quantify the finite size corrections. 



3 Numerical simulations 



The properties of the system are numerically studied via the Metropolis Monte 
Carlo algorithm [H] . We focus on the quenched scenario and reconstruct the 
average distribution of the main quantities of interest by averaging over several 
realizations of the graph of connections. The analytical solution of the annealed 
model is sketched in Appendix A. 
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The simulated system reproduces well the phase transition, the actual value of 
the temperature associated with symmetry breaking depending on the number 
of simulated spins. To test the scenario discussed in the previous Section, we 
first estimate (5c using the so-called Binder cumulant [H], defined as 

where (m^) and (m^) denote respectively the second and fourth moments of 
the magnetization. The Binder cumulant is computed for different values of the 
imposed temperature. These numerical experiments are repeated for distinct 
values of N , while keeping 7 fixed. The obtained profiles Un vs. T are reported 
in Fig. [T]for various values of 7. Notice that we have introduced a subscript N 
to recall that the plotted profiles are reconstructed from finite N calculations. 
The importance of the cumulant concept stems from the observation that 
curves corresponding to different all intersect at approximately the same 
temperature, which provides an estimate of the critical temperature Tc = l/(3c 
in the infinite limit. A direct inspection of the enclosed figures, suggests that 
for all values of 7 scanning the relevant interval (1,2), Tc = 2. This result is 
in agreement with the convergence to the mean-field limit irrespectively of 7. 
Clearly, the convergence to the mean-field solution is expected to be faster for 
larger values of 7. Indeed, Fig. [H panel a), which refers to the case 7 = 1.2, 
shows a less clear intersection in the interval of A^ covered by our investigations, 
when compared to similar plots depicted for larger values of 7. 

More interestingly, working at finite A^, one can monitor the magnetization 
and plot it as a function of the dilution parameter 7. Results are reported in 
Fig. [2] for two different choices of A'" (symbols): A tendency to asymptotically 
approach the mean-field reference (solid) line is clearly displayed, in agreement 
with the above scenario. Finite size corrections play however a crucial role, 
which deserves to be carefully addressed. 

Aiming at shedding light onto this issue, we use in the following two different 
analytical methods to estimate the finite A^ corrections. 



4 The replica method 

The first method relies on the replica trick. Starting with the calculations of 
Section [2], we now include the leading order finite A^ corrections. We start 
again from Eq. (fTT]) . which we recall here 

(^">.-Eexp(i:[|T.,4|;i;5-i^i;5]). (20) 
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Fig. 1. Binder cumulants Un as a function of temperature T, for 7 = 1.2, 1.4, 1.5, 1.7. 
The time averages are calculated over 4.10^ Monte Carlo sweeps. The error bars are 
estimated with the resampling technique using 10 sets of 1000 sweeps and calculating 
the associated variance. 



The largest neglected term in this expansion is 2J aJt^^ij ) including it would 
imply coupling three replicas and make the calculation technically more dif- 
ficult. However, this term becomes progressively more important as 7 ap- 
proaches 1 from above; this may possibly affect the accuracy of the prediction 
derived here, in this region of parameters. In the following we shall also neglect 
0{N~'^) terms in expression (120]) . which is certainly well motivated as long as 
7 < 2 (in a strict sense). 



Notice that the following relations apply 



and 



ijLj a a i^Lj a i a 



Nn 



(21) 
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Fig. 2. Magnetization versus 7 for (3 = 1. Symbols refer to the quenched numerical 
simulations for, respectively, N = 1000 (circles) and A'^ = 8000 (triangles). Dashed 
lines are guides for the eye. The solid line stands for the mean- field solution (|16p . 

E ^? = E E sts'^sls] = E ( E sts';s\s] - n) 

ii'j «7^J a,b i,j 

= E ( E S-S^ '-n'N = Y: Kb + '^'N, (22) 

a,b i a<b 

where use has been made of the definitions 



ma = Y: S:/N , (23) 

i 

mab = Y.StSl/N . (24) 

i 

Substituting ^ and ([22]) into Eq. ^ yields 
{Z-)j = E exp (pNY^ml - n(3 + ^-^N^^' ^ ml, + ^-^N^-'n) , (25) 

{SI} ^ a ^ a<b ^ ^ 

where the term scaling as v? has been dropped (recall that we shall be con- 
cerned with the limit n — > 0). 

The Hubbard-Stratonovich identity can now be invoked to rewrite the above 
exponentials involving and m^^ 



exp (pNrr?^ = J ^ J^^ dK exp (-^NXl + 2/3iVm,A,) , (26) 
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exp ^m^b = \ n / ^1<^b 



exp 



27r 



oo 



N^~"'qlb + l\f3^N^-"'mabqab ■ (27) 



Putting the various pieces together, the average rephcated partition function 
reads 



/oo ^ 



n dqab exp (-/3N ^ + 2/?iV ^ m,A, - ^^iV^"^ ^ + ^ip^N'^' ^ m^^g.^) 



^ C n n exp (-/3iV 5: - ^^iV^-^ ^ g^,) (28) 
^ exp (2/5iV 5^ m^Xa + 7!/?' AT^-t ^ m^.g^^) , 



where the normahzation C(7, /3, n, A^) = ^'^^^^^ ^ can be safely ignored 

in the forthcoming development. 

Let us focus now on the last sum appearing in Eq. (125]) . A straightforward 
manipulation leads to 



J2 exp (2(3N rriaXa + i\(3^N^~^ ^abQab 

{S°-} ^ o. a<b 



: E exp (2f3N E E f + 7!/5^iV-- E lab E ^) (29) 

{S^} ^ a i a<h i ^ 

E exp (2/3 E + 7!/9'Ari-^ E 9a6^'^5^) , (30) 

{5"} ^ « a<b ^ j 

where the index % can be removed, being replaced by the power A^. 
Let us now introduce the function as 



"^{Xa.qab) = E exp [2PNY^aS'' + j\p^N'-^Y<labS^S'') . (31) 
We now make the replica symmetric hypothesis which corresponds to setting 
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Xa = X and qab = q Va, b. After this ansatz, \1/ can be cast in the form 



{Sa} ^ a a 



The Hubbard- Stratonovich trick allows us to write 

^1 roo 



(32) 



exp (^iV^^ME^")'' 



(33) 



which leads to the following expression for 



^(A, q) = exp I I W — J^^ dx exp 



E exp ((2/3A + 2^m'^^x) ^ 5") 



■ exp 



71/32 A^i-Tgri\ /t! /•+00 

2 fei-oo ^^'"P 



7! 2 
• — X 

2 , 



7! 2 



1-7 



2 cosh(2/3A + 7!/? AT— ygx) 



Finally, we obtain 



where 



/oo 
dAdgexp[ArL„] 
-00 



(34) 

(35) 
(36) 



We want to keep the terms linear in n in Eq. (15^ , since they are the only ones 
to give a contribution in the limit n 0. A straightforward expansion in n 
yields, in the n — > limit. 



- In ^ ^ — -flp^N^^^q + In 2 + W — / exp ( -7!xV2 ) 
n 2 \ 2tt J-oo ^ ' 

In cosh (2/5A + ^\^N^^-^^I'^ ^x) dx . (37) 

The computation of the free energy is now reduced to finding the saddle points 
of the following function 



47! ' 27 

Incosh (2/5A + e^/^a;) dx 



(3J 
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Fig. 3. Magnetization versus 7 for (3 = 1. Symbols refer to numerical simulations 
for N = 1000 (circles). The dashed lines serves as a guide for the eye. The solid line 
stands for the replica based solution (f42l) . 



where we have introduced the small parameter 

e = 7!/5iV(i-^)/2 . (39) 

Expanding in e up to order and performing the Gaussian integrations, we 
get 

2 2 

0(A, q) = -(3\^ + In cosh 2/3A + - tanh^ 2(3\ + In 2 + oie^) . (40) 

47! 27! 

At order e*^, the conditions d\(f) = and dq(j) = yield 

Ao = tanh2/5Ao , go = tanh^ 2p\o . (41) 

As it should, the mean-field solution is recovered from these equations, see 
Eq. (fT6l) . We now write A = Aq + e^Xi, and we get from the condition d\(j) = 

_ 1 ^0(1 - ^0) /.^N 

- -^1- 2,3(1 -A§) ■ ('2' 

The dummy parameter A can be shown to correspond to the magnetization. 
Thus, one can compare the replica based prediction Eq. fll2]) (which takes into 
account the leading order finite corrections) to direct numerical simula- 
tions. The comparison is made in Fig. [3], where the magnetization is plotted 
as a function of 7. The global trend is captured by Eq. (l42l) : however the 
agreement deteriorates quickly for small values of 7. This may be due to the 
approximations involved in the calculation. 

As a final step, let us compute the leading finite A^ correction to the transition 
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temperature. We recall that its — oo value is /3c = 1/2. We compute the 
Hessian matrix H^{Q, 0) of in (A, q) = (0, 0), getting 



The critical temperature corresponds to a vanishing determinant for H^{{], 0). 
We find, at order e^, j3c = 1/2. In conclusion, at this level of approximation, 
there is no modification of the critical temperature due to finite effects. 

We end this Section with a comment. A given finite random graph with A^ sites 
and M links may be seen as a finite A^ realization of a dilute random graph as 
above, for some value of 7. It may also actually be seen as a finite A^ realization 
of a graph constructed according to the rule described in Eq. (jlj) and Eq. ([5]), 
with p = aNi/N ~ '^N'^~'^ . To each choice of 7 g]1, 2[ corresponds a value of 
a. There is then an infinite number of models to which our graph at hand may 
be compared. However, the result of a replica calculation at first subleading 
order in 1/A^, as performed above for a = 1, does not depend on the choice of 
7 and a. This freedom thus cannot be used to optimize our predictions for a 
finite A^ graph. 



5 An alternative approach: The cavity method 

The method used in the previous Section (expansion in powers of A^ coupled to 
a replica calculation) does not give very precise results for small to moderate 
values of 7. We now turn to an alternative theoretical approach to interpret the 
results of our simulations: a finite size graph with A^ sites, constructed with 
the rule ^ for a given 7, may be seen as a standard Erdos-Renyi random 
graph with parameter Xer = N^"'~^y'y\ [TT|3] . This makes possible the use of 
the powerful methods devised for finite connectivity random graphs, such as 
the cavity method. 

The solution of the Ising model on random graphs with arbitrary distributions 
of links is given in [T61IT7] ; we follow here the formulation given in [13] , which 
is there applied to the solution of an Ising spin glass on a Bethe lattices. Our 
case is much simpler, as we are studying a ferro magnet; a small complication 
is related to the probability distribution of the site connectivities. 

We briefly recall the main steps leading to the (replica symmetric) cavity equa- 
tions, following Ref. [13]. Consider the Hamiltonian ([1]) defined on a random 
graph. Figure m represents a node, denoted with 0, and its = 3 neighbours. 
We represent hj hi,i = 1, ...,k the total field acting on spin Si in the absence of 




(43) 
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Fig. 4. This figure shows an example of a tree-hke structure with k = 3. The cavity 
fields hi represent the total field acting on the spin Si, when the central spin Sq is 
removed. 

the central spin 5*0. The magnetization of the i-th spin reads = tanh(/3/ij). 
The basic ingredient of the cavity method is to assume that the fields hi are 
uncorrelated. 

Let us define the partition function of the spin 5*0 as follows 

k k 

Zs, = E exp {p'So E JSi + /?' E hiS^) , (44) 

So,Si,..,Sk i=l i=l 

where 

0' = 1^/3 . (45) 

One can now invoke the basic identity 

E exp (^(3'SoJSi + 0hiS^ = c{J, hi) exp(/5'u(J, hi)So) , (46) 

Si=±l 

where the two functions u{-,-) and c(-, ■) respectively read 

u{J,h) = — arctanh[tanh(/3 J) tanh(/5 h)] , (47) 

c(J./,.)^2 ^°-"'f'^>7f'f) . (48) 

cosh(p u{J, h)) 

and rewrite the partition function (jH]) as 

k 

Zs, = E niic( J, h) exp (/?' ^ m( J, /i,)^o) . (49) 

So,S\,..,Si^ 1=1 



In practice the magnetization on site is thus given by mo = (Sq) = tanh(/5'/io), 
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where 

fc 

ho = Y.u{J,h,) . (50) 

1=1 

The connectivity /c of a given site in the finite size random graph we study is 
a random variable with distribution 

Mk) = fc,(ir-"i-fc)! (l - N^-r-'N^^-'^' . (51) 

We want to compare our finite size random graph with the corresponding 
infinite size graph with the same Erdos-Renyi parameter Xer = N'^~^/'y\. 
This graph has a Poissonian connectivity distribution with paramter Xer- 
Thus, we use in the cavity calculations the following distribution, which is 
close to flSTl) . 

\ k 

Tr(k) = e-^^«^ . (52) 
k 

Eqs. (l50l) and ( l52l) allow one to write the following implicit relation for the 
probability density Q{h) of local fields 

k 

Q{h)=Y,TT{k) f Ul,[dh,Q{h,)]6{h - T.^{J^,k)). (53) 

k •' 1=1 

Eq. (l53l) can be solved using a population dynamics algorithm [13] to ac- 
cess an estimate of the local field distribution Q{h), and eventually compute 
the magnetization of the system. More concretely, one starts with arbitrarily 
chosen population of M fields and proceeds iteratively as follows. A random 
number k is picked up with probability 7i{k); a subset of k fields is randomly 
selected in the population, and used to compute the ho field, as prescribed 
by Eq. (I5UI) . Then one field is removed at random from the population, and 
replaced with the computed ho- Such a scheme defines a Markov chain on 
the space of the M fields which admits a stationary distribution. In the limit 
M — >■ cx) such a stationary distribution clearly satisfies the self-consistency 
relation 0531) . The magnetization m is hence straightforwardly recovered as 
m = X]jtanh(/? hi)/N. 

This method provides very accurate predictions, as shown in Fig. [5l much 
better than the finite expansion around the replica calculation discussed in 
Section HI 



6 Conclusions 

We have investigated the finite size corrections to the equilibrium magneti- 
zation of the Ising model defined on a diluted network. Varying the dilution 
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Y 

Fig. 5. Magnetization versus 7 for (3 = 1. Symbols refer to numerical simulations for 
= 1000 (circles). The solid line is the solution obtained using the cavity method. 
The agreement is much better than for the replica approach (see Figl3]). 

parameter 7, these networks interpolate between the fully connected network 
(7 = 2) and the opposite setting where the number of links scales linearly 
with system size (7 = 1). 

Systematic deviations with respect to the asymptotic mean field behavior are 
observed when a finite number of spins is considered, such discrepancies being 
more pronounced as 7 approaches its lower bound 7=1. This phenomenon is 
clearly displayed in the plot of the magnetization m versus 7. A replica based 
perturbative analysis is developed, whose predictions are compared with the 
outcome of the numerics. The dependence of m on 7 is qualitatively captured, 
but the quantitative match is not satisfying, especially as the dilution rate is 
increased. A cavity based calculation, inspired by the Mezard-Parisi technique 
[T3] is able to reproduce the data with an excellent degree of accuracy. 

Summing up, we have brought convincing evidences that finite size corrections 
do play an important role in presence of a diluted network and thus need to 
be carefully addressed. This is best done by using the cavity metod and the 
associated population dynamics algorithm. Further extension of the present 
analysis would include clarifying the reasons why the replica method turns 
out to be less accurate than the alternative cavity based approach. This may 
imply pushing further the replica calculation by acommodating for the so far 
neglected couphng among three independent replicas. 
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A Appendix A 



We shall give here the annealed solution of model ([T]). Averaging the parti- 
tion function Eq. ([2]) over the Jjj's, using the probability distribution [31 one 
straightforwardly obtains 

{Z)j =J2 n [1 - 4^^'' + ^N^-'e^'^'^] 
{s,} <id> ^- ^■ 

= ^ exp[^ln(l - ^N^-^ + l^N^-^e^^^'^')]. (A.l) 

{S,} i,j 

To proceed further we shall recall that 7 > 1, an observation which in turn 
allows us to expand in power of 1/A^ the above expression. The following 
expression is formally recovered 

{Z)J=T.^^P{^T.SA) , (A.2) 

{Si} i,j 

where the finite N temperature, (3^, reads 

In the above derivation we made use of the fact that (SiSj)"^ = 1 for m even 
and (SiSj)"^ = SiSj otherwise. In the limit for — ^ cxd, Eq. (1A.3I) implies 
Pn ^ P, which in turn implies 

(Z),= Eexp(^E^^^.) ^ (A.4) 

{Si} i,j 

for each value of the 7 parameter. The annealed solution is thus also equiv- 
alent to the fully coupled graph solution, in the N ^ 00 limit. For finite 
N, Eq. (1A.3P implies a modification of the temperature due to finite size ef- 
fects. According to (lA.Sp . we can imagine to replace the finite N system at 
temperature T = l/kB^ with its Curie- Weiss counterpart, provided a slightly 
smaller value of the temperature is allowed. This finding would in turn suggest 
that the finite graininess of the distribution drives an increase of the critical 
temperature. Consequently, one would expect to observe an inhomogeneous 
state at the mean-field transition temperature. This is at variance with what 
is found in our (quenched) simulations, where an opposite tendency is man- 
ifested. This discrepancy is also signaled by inspecting the magnetization as 
a function of 7, as outlined in Fig. lA.ll The numerics in Section [3] and the 
replica based analysis of Section [5l that both share the quenched viewpoint, 
display a similar trend (though the matching is not perfect as commented in 
the body of the paper). Conversely, the annealed prediction obtained from 
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Fig. A.l. Magnetization versus 7 for /? = 1. The solid line refers to the mean- field 
solution; the dashed line stands for the replica based calculation; the dotted line 
represents the annealed perturbative estimate. 



the mean-field magnetization associated to the finite N temperature (1A.3P 
returns a striking different behaviour. 



B Appendix B 



Let us here discuss the solution of the mean- field model. By using the Hubbard- 
Stratonovich transformation 



exp 



(ba"^^ = rfx exp + 2a6x) , (B.l) 

the partition function flA.4p (hereafter simply Z) can be cast in the form 



{S^} \ t / {Si} \ 



TX 



/■"J rfxexp {-fiNx') |: exp (2/3A^E 
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where a = (J2i SiY and b = PN and eventually 



Z 




oo 



dxexp (-iV(/3x^ - ln(2cosh(2/3x)))) 



(B.3) 



The free energy function results in 



-(3F= lim — InZ, 



N — 




N 71 J-O. 



■OO 




ln(2 cosh(2/3 
(B.4) 



and using the saddle point approximation 



pF = max[-/3x^ + ln(2 cosh(2/3x))]. 



(B.5) 



References 

[I] A. Barrat, M. Barthelemy and A. Vespignani, Dynamical processes on complex 
networks, Cambridge University Press (2008). 

[2] S. N. Dorogovtsev and J. F. F. Mendes, Evolution of networks, Oxford 
University Press (2003). 

[3] S. N. Dorogovtsev, A. V. Goltsev and J. F. F. Mendes, Rev. Mod. Phys, 80, 
1275 (2008) 

[4] A. K. Hartmann and M. Weigt, Phase transitions in combinatorial optimization 
problems, Wiley (2005). 

[5] L. De Sanctis and F. Guerra, J. Stat. Phys., 132, 759 (2008). 

[6] A. Barra, J. Stat. Phys., 132, 787 (2008). 

[7] A. Barra and P. Contucci, larXiv:0812.T568] (preprint 2008). 

[8] E. Aghari, A. Barra and F. Camboni, arXiv:08044503 (preprint 2008). 

[9] A. Dembo and A. Montanari, arXiv:08044726 (preprint 2008). 

[10] A. Bovier and V. Gayrard, J. Stat. Phys., 72, 643 (1993). 

[II] P. Erdos and A. Renyi, Pubhcationes Mathematicae 6 290-297 (1959). 

[12] M. Mezard, G. Parisi and M. A. Virasoro, Spin glass theory and beyond, 
Singapore, World Scientific (1987). 

[13] M. Mezard and G. Parisi, Eur. Phys. J. B, 20, 217 (2001). 

[14] N. Metropolis et al., J. Chem Phys., 21, 1087 (1953). 

[15] K. Binder, Z. Physik B 43, 119 (1981); Phys. Rev. Lett. 47, 693 (1981). 



20 



[16] S. N. Dorogovtsev, A. V. Goltsev and J. F. F. Mendes, Phys. Rev. E, 66, 016104 
(2002). 

[17] M. Leone, A. Vazquez, A. Vespignani and R. Zecchina, Eur. Phys. J. B, 28, 
191 (2002). 



21 



